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We present Feynman type diagrams for calculating the n-point function of the primordial curva- 
ture perturbation in terms of scalar field perturbations during inflation. The diagrams can be used 
to evaluate the corresponding terms in the n-point function at tree level or any required loop level. 
Rules are presented for drawing the diagrams and writing down the corresponding terms in real 
space and Fourier space. We show that vertices can be renormalised to automatically account for 
diagrams with dressed vertices. We apply these rules to calculate the primordial power spectrum 
up to two loops, the bispectrum including loop corrections, and the trispectrum. 

PACS numbers: 98.80.Cq YITP-07-30 



INTRODUCTION 



There is currently a great deal of interest in the statistical properties of primordial perturbations from inflation, 
, because measurements of any non-Gaussianity will improve by about an order of magnitude over the next few years, 
^ ' for example with the Atacama Cosmology Telescope [l| and Planck ^ . This will provide a key way to discriminate 
, between the many models of inflation. Although single field models of slow-roll inflation typically generate a small 
level of non-Gaussianity_[l, 01, there may be an observable level generated in multiple field inflation [EL for example 
' the curvaton scenario K|, or some single field models such as Dirac-Born- Infold inflation (DBI) [3]. Currently 
^ ' observations of the CMB have concentrated on constraining the 3-point function (bispectrum) [3, j but the 
\^ \ 4-point function (trispectrum) has also been considered [ll|, [Tsj , and in principle higher order n-point functions 
0^ ■ may also be observable if they are sufficiently large. In some models the first signal of non-Gaussianity may come 
Op ' through the trispectrum, for example some special cases of the curvaton scenario [14]. Furthermore higher order 
statistics in principle carry more information and therefore could distinguish between different sources of non- 
10 ' Gaussianity, for example in many cases the bispectrum can be paremetrised by a single non-linearity parameter, 
! Jnl, while the trispectrum depends on two wavenumber independent parameters, tjvl and qnl [T3 |. 
■ We will calculate the primordial curvature perturbation on uniform density hypersurfaces, C, on large scales 
employing the 5N formalism [H, [l^, [13, [ll] . We use the separate universe approach [H, [l^, [2l| . This considers 
each super-Hubble scale patch to be evolving like a separate Friedmann-Robertson- Walker universe which are locally 
homogeneous. By patching these regions together we can track the evolution of the perturbations on large scales just 
by using background quantities which greatly simplifies our calculations. We also need to know the perturbations 



d of the scalar fields at Hubble-exit, but in many cases these are extremely close to Gaussian in which case their 
statistical properties are given purely in terms of the power spectrum of the fields. 
The number of e-foldings, iV, given by 

N = j^^ "° H{t)dt , (1) 

is evaluated from an initial flat hypersurface to a final uniform-density hypersurface. The perturbation in the number 
of e-foldings, 5N , is the difference between the curvature perturbations on the initial and final hypersurfaces. We 
wish to calculate primordial perturbations, hence we pick a final uniform density hypersurface to be at a fixed time 
during the standard radiation dominated era, for example during primordial nucleosynthesis. The initial time is 
arbitrary provided it is after the Hubble exit time of all relevant scales. It is often convenient to pick this time to be 
shortly after Hubble-exit. We introduce the method for calculating the primordial curvature perturbation in terms 
of the 5N formalism in more detail below. 

Here we present a Feynman diagram type approach to calculating the primordial non-Gaussianity, where every 
term in the n-point function of Q has a diagrammatic representation. Similar diagrams were introduced in the 
context of non-linear perturbation theory of large scale structure [l^l and further developed in [2^ (see also f23|). 
It is possible to read off the corresponding mathematical term in the n-point function for every diagram and it is 
possible to draw all the diagrams in a systematic way so that all terms are included. This approach extends to loop 
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corrections of any required order. Although the loop corrections are expected to be small in general compared to the 
tree level terms there is no proof of this, except for a particular loop correction to the bispectrum in some specific 
models of inflation [2^ . We will see that in general the loop corrections depend on a large scale cut off; de pen ding 
on this cut off the loop correction may contribute significantly to observable quantities, see for example [251 [26j. 
In fact there is a specific model of inflation where the one loop term gives the dominant and potentially observable 
contribution to the bispectrum '27\. A diagrammatic method to calculate field fluctuations during inflation was 
presented in psj . 

We first introduce the primordial curvature perturbation. In the next two sections of this paper we assume the 
initial field perturbations are Gaussian, which in many cases is a good approximation to make, for example slow-roll 
inflation. In Sec. IIII Al the rules for drawing diagrams in real space are given, along with a discussion of why we 
want to calculate connected rather than disconnected diagrams. In Sec. IIII 5] the equivalent rules in Fourier space 
are given. An application of the Fourier space rules for the power spectrum including second order loop corrections 
is given in Sec. IIII Cl In Sec. lIVI we discuss a way to renormalise the vertices of the diagrams such that all diagrams 
with dressed vertices can be absorbed, with a proof given in the appendix. In Sec. |V]the assumption of Gaussian 
initial fields is dropped and the complete rules presented and applied to several examples. Finally we conclude in 
Sec. ED 



II. THE PRIMORDIAL CURVATURE PERTURBATION 



We can write the primordial curvature perturbation, C, in terms of derivatives of N with respect to the fields, 
Na = dN/d(j)"^ and the field perturbation, ip^, at the initial time tini- By neglecting any additional dependence on 
(j)^ at Hubble-exit we are assuming the field perturbations on large scales are overdamped and thus we can neglect 
the decaying mode. This is valid during slow roll, but also more generally holds for light fields during inflation. The 
assumption that iV is a function of field values only is also valid in DBI inflation We use the notation that 
the fields are cp^ — cf)^ + ip-^, where the perturbed field ip"^ = Scj)"^ satisfies (ip^) = 0. The curvature perturbation 
is given by 

C ^ SN- {SN) = Na^^ + \nab iv^p'' - (^^^^)) + I^Nabc {^^v''^'' - i^^v''^'')) + • ■ • , (2) 
where SN = Na^^ + ^Nab^^P^ + ^NabcP^V^P^ + ■■■ ■ (3) 

The dummy index A labels the light scalar fields relevant during inflation, and summation convention is used 
throughout this paper. It is more convenient to work with the C defined above since it satisfies (C) = 0, even though 
often C — SN is used in the literature. This equation can hence be used to calculate the primordial n-point function 
of C, although the result will depend on the n-point function of the fields at Hubble exit. This has so far been 
explicitly calculated for the 2-, 3- and 4-point functions, [4, 30., ,31]. 
The connected 2-, 3- and 4-point functions of the fields are defined by 

«0 - C^^{k){2nfS'{k^+k2), (4) 
{<<P^,) = i?^^^(fci,fc2,A:3)(27rf<53(ki+k2+k3), (5) 

(y^k.^k.^^Oc = T^^^^(ki,k2,k3,k4)(2^f<53(ki+k2+k3+k4). (6) 

Note that only the 4-point function (and higher) depend on the direction of the k vectors, the 2- and 3-point 
functions of the fields just depend on the magnitude of the vectors, ki = |ki|. Homogeneity of the random fields 
implies that the sum of the k vectors must be zero, while isotropy implies that the n-point functions are invariant 
under reorientations of this closed configuration. 

At lowest order in slow roll different fields are uncorrelated at Hubble exit and all light fields have the same power 
spectrum [H, H^l , so 

C^B{k) = S^^P{k), (7) 
where S^^ is the Kronecker delta-function, and the variance per logarithmic interval in fc-space is given by 

where the Hubble parameter H is evaluated at Hubble-exit, k ~ (ai?)*. At zeroth order in slow-roll parameters, 
V is independent of wavenumber, i.e. we have a scale invariant spectrum for the field fluctuations. The bispectrum 
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satisfies B^^^ = 0{e^/'^P^) where e is a slow-roll parameter, so it is zero at lowest order in slow roll. However 
the trispectrum is not zero at lowest order in slow roll [31[, T^^^^ = 0{P^)^ but is suppressed by 0{P^) compared 
to the power spectrum. For a discussion of the slow-roll order of higher order n-point functions of the fields see [3^] . 

III. DIAGRAMS FOR GAUSSIAN FIELD PERTURBATIONS 

Initially we will assume that all of the fields have a purely Gaussian distribution at an initial time, e.g. shortly 
after Hubble exit. Then every odd n-point function of the fields is zero and every even n-point function of the fields 
can be reduced to a product of 2-point functions. This assumption is not required and the extension to non-Gaussian 
field perturbations will be presented in Sec.|Vl However the approximation of initially Gaussian field perturbations 
is frequently made in the literature and is often a good approximation, for example [s^ shows that the 3-point 
function of the fields adds an unobservably small contribution to the bispectrum around Hubble exit assuming slow 
roll inflation with a standard kinetic term of the scalar fields, while [31(1 show that the 4-point function of the fields 
adds an unobservably small contribution to the trispectrum around Hubble exit under the same assumptions. 

Working to zeroth order in slow roll implies that ([7]) holds and that the 3-point function of the fields is zero. 
So although working to lowest order in slow roll is not the same as working with Gaussian field perturbations, in 
practice the conditions are related. Hence when assuming that the field perturbations have a Gaussian distribution 
at Hubble exit we will also work at zeroth order in slow roll. It is likely that any observable amount of non- 
Gaussianity can be calculated to sufficiently high accuracy when making these two approximations, assuming a 
standard kinetic term. 



A. Real space diagrams 

Assuming that the fields arc uncorrelated Gaussian variables, with identical distributions, as discussed above, the 
two point function of the fields is given by 

, ARx_{ 5'^^G{\x^-Xj\) fori^j, , . 

('■^x.'^.,) - I ^AS^^2^ fori = j. 

All of the higher n-point functions of the fields can be written in terms of G and the variance of the fields 
The two point function G only depends on \xi — Xj \ because we are assuming the background is homogenous and 
isotropic. 

The diagrammatic rules for the connected n-point function, {C,xiC,x2 ' ' ■Ca;„)c, are: 

1. Draw n points representing the n spatial points xi, - ■ ■ ,a;„ and connect them with r propagators (dashed 
lines) which attach two of the positions Xi (it can be two distinct positions or it can have both ends attached 
to the same position) . The diagram should be connected in order to calculate the connected n point function, 
see Sec. IIII A II We require r > n — 1 in order to draw a connected diagram and diagrams with r = n — 1 are 
tree level, while those with r > n — 1 include loop corrections. 

2. Label each end of each propagator with the field indices A, B, . . . C. 

3. Assign a factor Nab---c to each spatial point, x^, where the number of derivatives of N is the number of 
propagators attached to that Xi. 

4. Assign a factor of 5^^G{\xi — Xj\) to each propagator, where AB are the appropriate field indices attached 
to the propagator and Xi, Xj are the positions at either end of the propagator. If both ends of the propagator 
are attached to the same Xi instead assign a factor 5^^{(p^). 

5. Divide by the appropriate numerical factor. Whenever / propagators attach the same Xi and Xj at both ends 
this gives a factor of /!. Whenever I propagators dress an Xi this gives an additional factor of 2' (as well as the 
factor of l\ due to the previous rule). A propagator dresses an Xi if both ends of the propagator are attached 
to the same Xi. 

6. Add all permutations of the diagrams which is all of the distinct ways to relabel the spatial points. The number 
of permutations depends on the symmetries of the diagram, a diagram with complete symmetry between all 
of the spatial points has only one term, while for a diagram with no symmetries between the spatial points 
there are n! permutations. 
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FIG. 1: The 2-point function at tree (a) and one loop level (b), and the 3- and 4-point functions at tree level, (c) and (d) 
respectively. All figures in this paper were drawn using JaxoDraw [s^. 



In Fig. [T] the diagrams for the 2-point function at tree level and one loop level, plus the tree level terms for the 
3- and 4-point functions are shown. Note that for the 4-point function there are two tree level terms. 
The terms corresponding to the diagrams are given below, in the same order as the diagrams, 

(GiC.) = {NaN^ + NaN^''{^''))G{\xi~X2\) + \nabN^''G{\xi-X2\)\ (10) 

iUCx^U) = NaNeN^"" (G(|xi - X2\)G{\xi -x^\) + 2 perms) , (11) 
{Cx,Cx,Cx,Cx,)c = NaNbNcN^''^ (G(|xi - a;2|)G(|:Ei - X3|)G(|a;i - X4I) + 3 perms) 

+NaNbN^N^^ {G{\xi - a;2|)G(|x2 - X3\)G{\x3 - Xi\) + 11 perms) . (12) 

Although the second diagram for the two point function. Fig. [1] has a dressed vertex with an associated numerical 
factor of 2, there is another permutation of this diagram with the X2 dressed which gives an equal contribution. The 
third diagram has two propagators attaching xi and X2 and there are no other permutations of this diagram so there 
is a numerical factor of 2, as shown in the third term of pH)) . The diagrams for the 3- and 4-point functions are all 
at tree level so they all have a numerical factor of 1 . The first diagram for the 4-point function has 4 permutations 
since we have 4 choices of which of the xi, • • ■ ,2:4 has 3 propagators attached to it as given in the first term of HI]), 
the second diagram of the 4-point function has 12 permutations because there are 4x3 choices of which 2 of the 
xi , • ■ • ,0:4 should have 2 propagators attached to them 



1. Disconnected diagrams 

In general any n-point function has a connected and a disconnected part. Throughout this paper we will only 
consider the connected n-point function because the disconnected contribution contains no new information over 
the lower n-point functions that it is split into. There is no disconnected contribution to the 2- or 3-point function 
of the fields or zeta, because we assume = (C) = 0. The disconnected contribution is only non-zero if a subset 
of the k vectors sums to zero. 

For example the disconnected part of the primordial 4-point function, (Cj^^ CkaCk^)' ^^"^ product of two 
2-point functions such as (^^^ Cjfj ) (Ckj Ck^ ) ^^-^ in total there are three such terms. This is non-zero even if C, is 
purely Gaussian, but it contains no new information compared to the primordial 2-point function. This term only 
contributes if -I- kj =0 for some i,j = 1, • • • ,4, so provided the sum of any two k vectors is not zero any 
contribution is from the connected part. The connected part of the primordial 4-point function is only non-zero 
if C, is not purely Gaussian and it contains information that is not observable from the primordial 2- or 3-point 
functions. 



B. Fourier space diagrams 

We present here the rules for drawing the connected diagrams (see Sec. IIII A 1[) of the n-point function of C, at 
r-th order (i.e. 0{P'^)), for r > n — 1. The tree level terms correspond to r = n — 1. 
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FIG. 2: The terms that are associated with every vertex and propagator. 

1. Draw all distinct connected diagrams with n-external (solid) lines and r (dashed) propagators. Every vertex 
must consist of 1 external line and at least 1 propagator. 

2. Label the external legs with incoming momenta for i — 1, • ■ • ,n and label the propagators with internal 
momenta for i = 1, ■ • ■ , r. Label each end of each propagator with a field index A, B, ■ ■ ■ ,C. 

3. Assign a factor Af^B - c(27r)'^(5'^(ki — qi — • • • — qp) to each vertex. The number of derivatives of N corre- 
sponds to the number of propagators attached to each vertex. We use the convention that incoming momentum 
is positive. The d function ensures momentum is conserved at each vertex. See Fig. [S] 

4. Assign a factor 5^^ P{q) to each propagator, where AB are the appropriate field indices that the propagator 
is labelled with at either vertex and q is the momentum attached to the propagator. See Fig. [2l 

5. Integrate over the propagator momenta, -j^^p^ J (pQi- The first ri — 1 integrations are trivial because of the S 
functions but any further integrations (in the case of a diagram with loop corrections) cannot in general be 
evaluated analytically. 

6. Divide by the appropriate numerical factor. Whenever I propagators attach the same vertices at both ends 
this gives a factor of II. Whenever / propagators dress a vertex this gives an additional factor of 2' (as well 
as the factor of 11 due to the previous rule). A propagator dresses a vertex if both ends of the propagator are 
attached to the same vertex. 

7. Add all permutations of the diagrams which is all of the distinct ways to relabel the kj attached to the 
external lines. The number of permutations depends on the symmetries of the diagram, a diagram with 
complete symmetry between all of the external lines has only one term, while for a diagram with no symmetries 
between the external lines there are n! permutations. 
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1. Comments on the diagrams and their rules 



Going between the unlabelled diagrams in Fourier or real space is straightforward, to go from Fourier space to 
real space one simply removes the external lines and places an Xi at each vertex. To go from real space to Fourier 
space one attaches an external line in place of every Xi. 

In Fourier space rule 1 for the diagram is to draw every possible diagram at the appropriate order. Rule 2 gives 
the rules for labelling the diagrams. Rules 3-5 give the associated mathematical expression for each diagram, up 
to a numerical factor given in rule 6. A derivation of the numerical factors is given in Sec. IIVI and App. lA II The 
numerical factor for all tree level diagrams is 1. Finally rule 7 tells one to include every permutation of the labelling 
of the ki on the external legs that are distinct. 
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FIG. 3: The tree level (a) and 1 loop correction (b) for the power spectrum. 
' ) ' ) ' ) 



I ■ ■ ^ 

FIG. 4: The two loop terms for the power spectrum with a Gaussian initial field. 

C. Power spectrum 

The primordial power spectrum in Fourier space is defined by 

(Ck,Ck.)^Pc(fc)(27r)3<5^(ki+k2). (13) 

The tree level and one loop correction terms to the 2-point function is diagrammatically given by Fig. [31 After 
carrying out the trivial integration over the propagator momentum the tree level term is given by 

P^t>^<=<=(fc) = NAN'^P{k) . (14) 

The one loop terms are given by 

Pi -ji^Jd'q QiVABiV^^P(<z)P(|ki - q|) + NANi^Pik)P{q)^ , (15) 

where we have already carried out one integration. Both of the diagrams associated with these terms have a 
numerical factor of 2, however there are also 2 permutations of the second term because it has a non symmetric 
diagram. Because of momentum conservation, ki + k2 = 0, it follows that ki — k2. After enforcing k = ki = k2 
the two permutations of the second diagram give an equal contribution which cancels the numerical factor of 2. 

Going to 2 loops, terms of order there are four terms. Diagrammatically the terms are given in Fig. 2] which 
corresponds to 

= y'dVrf'g2(^A^;^lcA^^P(9i)P(<Z2)P(fc) + iA^iBA^<f^^(fc)^('zi)^(92) 

+ ^NABN^''^P{qi)P{q2)P{\ki - q^l) + ^NabcN^""^ P{qi)P{\ci2 ~ qi|)P(|ki - qsl)^ . (16) 

Only the fourth diagram of Fig. |4] does not have a dressed vertex. In fact at every level there is only one diagram 
without dressed vertices, as has already been shown explicitly at tree and one loop level in Fig. O 

D. Evaluation of the 1 loop diagrams 

In the case where a single field (j) generates the primordial curvature perturbation, so N — N{(j)), we have 

Pl'^'^k) ^ {N'fP{k) , (17) 

where TV' = dN/ dcj). In the case of a scale invariant spectrum we can also evaluate the 1 loop integrals if we apply 
a large scale cut off L. For a discussion of how to evaluate the integral over the loop momenta and the dependence 
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of this term on the cut off L see for example ^f\. We also apply a small scale cut off k. The result is 

pI i°°P(fc) ^ P{k)V {{N" f + N"'N') log(fcL) , (18) 

where the variance V was defined by Eq. ([5]). In analogy with Eq. ([5]) we define Vc^ = P(;(fc)fc'^/(27r^) which from 
observations of the CMB is of the order 10"^" [s^. We hence have 

We can write the derivatives of N in terms of the (in principle) observable non-linearity parameters /jvl and qnl 
which in this case can be defined by [13, [3§| 

C = Ci + IInlC! + ^9nlC! + ■■■, (20) 

where is Gaussian because it is directly proportional to the initial Gaussian field perturbation, (pi, and the 
dimensionless non-linearity parameters, f^L and gNL: are given by 

25 N'" 

- 54(W ^ ^ 



For the extension of these two formula to the case of multiple field inflation see for example . Substituting (|2T 
and (1221) into IHl) we find 



ptree(fc) 'C 1^25-'^^ 25' 



= ^fNL + 7T75iVL log(A;L) . (23) 



If 75*"^°° 10 and we take the large scale cut off to be comparable to the present day Hubble scale then we can 
take log(A:L) = 0(1) ^,[131, and the observational bound on the bispectrum, I/atlI ^ 100 [SS*! bounds the first 
term of the ratio above to be less than 10~^. There is not yet any equivalent observational bound on g^i but it 
would have to be extremely large in order to give a significant contribution to Eq. (|23p . It therefore seems that the 
only way for the 1 loop contribution to be significant compared to the tree level term of the power spectrum is to 
take an exponentially large cut off L. For a discussion of this possibility see for example [25, 26]. 



IV. RENORMALISATION 



There is a way to reduce the number of diagrams significantly, by renormalising the vertices such that the diagrams 
with dressed vertices are automatically accounted for. We do this by renormalising the factors Nab---c, that are 
attached to each vertex. 

The derivative of the number of e-foldings for the given background (/)o is Nab---c = ^AS - cUoi ^-nd we can 
relate this to the number of e-foldings at a general point x by 

N EE iV(0(x)) = N\^, + Nav^ + \nabV^v'' + j^Nabc^^V^'v^ + • ■ ' ■ (24) 

The expectation of Nab---c is given by 

{Nab-c) = Nab-c + \n^b-cd{^^) + \KE-CDE{^^f + ■ ■ ■ , (25) 

where we have used Wick's theorem to decompose the expectation of powers of according to 

(^2") = (2n-l)!!(^2^", {if^--')^Q. (26) 

The double factorial is defined by (2n — 1)!! — {2n — l){2n — 3) • • ■ 5 • 3 • 1. The coefficient of the general term in ((25)) . 
((^^)"' is given by ([M)) and (|26p and is (2n — l)!!/(2n)! = l/(2"n!). This explains the numerical factor associated 
with every dressed vertex, see rule 6, Sec. IIIIBl 
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FIG. 5: Graphical representation of {Na}- 



In Fourier space the ensemble average of Nab---c is 

(Nab-c) - Nab-c + l^AB-CDj^ J d^QPil) + I^aE-cdej^ J d^qid^q2P{qi)P{\ci2 - qil) + • • • ,(27) 

which in real space is given by Eq. (|25p . 

In Fig.Elthe graphical representation of (Na) — Na + \N^g{ip'^) + \N^^(j{ip'^)'^ + • ■ • is given up to 2 loops. So 
far the derivatives of N had been calculated at 0o, however if we instead calculate the expectation of the derivative 
of at a general point, and attach this factor {Nab---c) to each vertex (which we call a renormaliscd vertex) then 
all diagrams with dressed vertices are automatically included in the expansion (j25p and don't need to be drawn. 
We provide a proof of this in appendix [XI 



A. Renormalised real space diagrams 



The new diagrammatic rules for the connected n-point function with renormalised vertices are: 

1. Draw n points representing the n spatial points xi, ■ ■ ■ ,x„ and connect them with r propagators (dashed 
lines) which attach two distinct positions Xi and Xj. The vertices are drawn with a solid dot to show that 
they are renormalised, see Fig. [5] for examples of such diagrams. The diagram should be connected in order to 
calculate the connected n point function, see Sec. IIII A II We require r > n — 1 in order to draw a connected 
diagram and diagrams with r = n — 1 are tree level, while those with r > n — \ include loop corrections. 

2. Label each end of each propagator with the field indices A^B , . . . C . 

3. Assign a factor {Nab---c) to each spatial point, x^, where the number of derivatives of N is the number of 
propagators attached to that Xi. 

4. Assign a factor of 5^^G{\xi — Xj\) to each propagator, where AB are the appropriate field indices attached 
to the propagator and Xi, Xj are the positions at either end of the propagator. 

5. Divide by the appropriate numerical factor. Whenever / propagators attach the same Xi and Xj at both ends 
this gives a factor of /!. 

6. Add all permutations of the diagrams which is all of the distinct ways to relabel the spatial points. The number 
of permutations depends on the symmetries of the diagram, a diagram with complete symmetry between all 
of the spatial points has only one term, while for a diagram with no symmetries between the spatial points 
there are n! permutations. 

In Fig. [5] the diagrams for the 2-point function up to two loop level is presented. 

The terms corresponding to the diagrams are given below, in the same order as the diagrams, 

{UU) = {NA){N^)G{\x^-X2\) + \{NAB){N^^)G{\x,~X2\f 

+ ^(iV^Bc)(iV^^'^)G(|a;i - X2\f + ■■■ . (28) 

Note that with the renormalised diagrams there is only one diagram at every order in loops. This result is also 
derived to all orders in App. |^ 
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(a) (b) (c) 

^Xg #X2 

FIG. 6: The 2-point function at tree, 1 loop and 2 loop levels ((a), (b) and (c) respectively), with renormalised propagators, 
(a) (b) ^ (c) _ 

— — — #^ y — — — 

FIG. 7: The power spectrum at tree, (a), 1, (b), and 2, (c), loop levels with renormalised propagators. 

B. Renormalised Fourier space diagrams 

Rules for drawing the connected diagrams (see Sec. lIII AT]) of the rt-point function of C at r-th order (i.e. 0{V^')), 
for r > n — 1. The tree level terms correspond to r = n — 1. 

1. Draw all distinct connected diagrams with n external (solid) lines and r (dashed) propagators. Every renor- 
malised vertex must consist of 1 external line and at least 1 propagator, and drawn with a solid dot to show 
that the vertex is renormalised. The propagators cannot have both ends attached to the same vertex. 

2. Label the external legs with incoming momenta for i = 1, ■ • ■ ,n and label the propagators with internal 
momenta for i = 1, • • • , r. Label each end of each propagator with a field index A, B, - ■ ■ ,C. 

3. Assign a factor {NAB---c)S^O^i Hi ■ ■ ■ Qp) to each vertex. The number of derivatives of N corresponds 
to the number of propagators attached to each vertex. We use the convention that incoming momentum is 
positive. The S function ensures momentum is conserved at each vertex. 

4. Assign a factor S^^P{q) to each propagator, where AB are the appropriate field indices that the propagator 
is labelled with at either vertex and q is the momentum attached to the propagator. 

5. Integrate over the propagator momenta, I ^-^li- The first n ~ 1 integrations are trivial because of the 5 
functions but any further integrations (in the case of a diagram with loop corrections) cannot in general be 
evaluated analytically. 

6. Divide by the appropriate numerical factor. Whenever I propagators attach the same vertices at both ends 
this gives a factor of /!. 

7. Add all permutations of the diagrams which is all of the distinct ways to relabel the attached to the 
external lines. The number of permutations depends on the symmetries of the diagram, a diagram with 
complete symmetry between all of the external lines has only one term, while for a diagram with no symmetries 
between the external lines there are n\ permutations. 

C. Power spectrum 

We again calculate the power spectrum up to 2 loop level, this time with renormalised vertices. There is only 
one diagram at every loop level, we present the three lowest order diagrams in Fig. [71 which correspond to 

1 1 

-}F{k) + 

1 1 



3! (27r)e 



J d^qid^q2{NABc){N^'''')P{qi)P{h2 - qi|)P(|q2 - ki|) . (29) 
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(b) 



FIG. 8: The bispectrum at tree, (a), and 1 loop, (b), levels, with renormalised propagators. 



D. Bispectrum 

We calculate the bispectrum including the leading order loop corrections. The bispectrum is defined by 

(CkiCk.Cks) =Sc(fci:^2,A:3)(2^f<53(ki+k2+k3). (30) 

At tree level there is one term and at 1 loop level there are two terms, which are shown diagrammatically by Fig. [8] 
and these give the contribution 

= {Nab){N^){N^) {P{ki)Pik2) + 2 perms) , (31) 

Sc'°°' = Jd'q(^{NAB){N^){N^^)P{q)P{\k^-Cl\)P{\k,+Cl\) 

+ l{NABc){N^){N'''')iPih)Piq)P{\k2-q\) + 5perms)] . (32) 



V. EXTENSION TO NON-GAUSSIAN FIELDS IN FOURIER SPACE 



If we don't assume that the scalar fields are Gaussian at the initial time (e.g. Hubble exit) we have to take account 
of the higher n-point functions of the fields. This is likely to be important in models with a non-standard kinetic 
term where the non-Gaussianity of the fields may be large at Hubble exit [13], for example DBI inflation inspired 
from brane world models, 0, [2§|. The connected 2-, 3- and 4-point function of the fields are defined in ([4])-([6|). 

Our analysis relies on the 5A^-formalism which assumes that is a function of the field values (f>^ but not the 
field derivatives independently as well. However we can still go beyond leading order in a slow-roll expansion and 
derive expressions which are valid to arbitrary order in slow roll, for example the different fields are in general 
correlated at Hubble exit, at first order in slow roll, see [s^ [ssj. 



1. Extended rules for the diagrams with non- Gaussian fields 



The diagrams are more complicated since while the 2-point function is still drawn as a propagator the higher 
n-point functions are drawn by an n-point vertex of dashed lines. 
The rules for the diagrams are now: 
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1. Draw all distinct connected diagrams with n-external (solid) lines and the appropriate number of internal 
correlators, which can be propagators (2-point functions) or higher n-point correlators of the fields. The order 
of the connected n-point function of the field perturbations is n — 1, i.e. {f^)c ~ 0{'P"~^). There cannot be 
more than one external line attached to each vertex. Every dashed line must be attached to an external line 
on at least one end, if it is attached to an external line on both ends it is a propagator (a correlator of order 
2). 

2. Label the external lines with incoming momenta for i — 1, ■ • • ,n and label each dashed line with internal 
momenta for i = 1, • ■ • , r. Label each dashed line at a vertex containing an external line with a field index 
A,B,---C. 

3. Assign a factor Nab - c{'^'^)^S^O^ — qi — ■ • • — qp) to each vertex which includes an external line, where the 
number of derivatives of N corresponds to the number of dashed lines attached to that vertex. The 6 function 
ensures momentum is conserved at each vertex. See Fig. [9] 

4. Assign a factor C^^{q) to each propagator, where AB are the appropriate field indices that the propagator 
is labelled with at either vertex and q is the momentum attached to the propagator. See Fig. [HI 

5. Assign a factor (27r)'^5^(qi + q2 + <l3)B^^'^'{qi,q2,q3) to each vertex of three dashed lines and no external 
lines (which corresponds to a 3-point function of the fields), where A, B and C are the field indices attached 
to the dashed lines where they meet a vertex with an external line. See Fig. [9] again. Similar rules hold for 
vertices of dashed lines at higher order. 

6. Integrate over all internal momenta, I '^^li- ^'^^ ^ ^^^^ diagram all of the integrations are trivial, but in 
the case of a diagram with loop corrections there will be some integrals which in general cannot be evaluated 
analytically. 

7. Divide by the appropriate numerical factor. Whenever I propagators or higher order correlators of the fields 
attach the same vertices (and with the same number of lines to each vertex) at both ends this gives a factor 
of If u legs from a higher order correlator are attached to the same vertex this gives an extra factor of u\. 
It follows that if I correlators, each of order p, all dress the same vertex then this gives a numerical factor 
of /!(p!)'. See App. lA 21 for a proof of these rules and see Fig. [13] for two example of dressed vertices, the 
associated terms with the numerical factors are given in (jA21[) . (|A22[) . 

8. Add all permutations of the diagrams which is all of the distinct ways to relabel the attached to the 
external lines. The number of permutations depends on the symmetries of the diagram, a diagram with 
complete symmetry between all of the external lines has only one term, while for a diagram with no symmetries 
between the external lines there are n\ permutations. 



A. Power spectrum 

We give again the power spectrum defined by (|13p . but this time for non-Gaussian initial fields and valid to all 
orders in slow roll. There is one extra diagram at one loop level, and six extra diagrams at two loop level. Because 
of momentum conservation, ki + k2 = 0, it follows that ki = k2 = k. Therefore even the non-symmetric diagrams 
which would have two permutations only correspond to a single term after enforcing this symmetry. 

The tree level term is given by 

pt™(fc) = NANeC^^ik) , (33) 

notice that at leading order in slow roll using ([7]) this reduces to ([H]). 
The one loop terms are given by 

= J d'q(^lNABNcDC^''{q)C^'' -dl) + NANBCDC^''ik)C^''iq) 

+NANBcB^''^{k, q, |ki - q|)^ , (34) 

where the extra one loop term which does not appear for Gaussian fields is shown in Fig. 1101 ^^'^ the first two terms 
were shown diagrammatically in Fig. 3] Going to 2 loops, terms of order there are ten terms, only the last four 
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/ e ''2 



(2ji)353(k-<|^-q2^l3)N 



3''^ ABC 



(A) 

'''' 



(2x)383(q^+q2+q3)BABC(q^,q^,q^) 



FIG. 9: The terms that are associated with every vertex, propagator and higher order correlator. 



\ > 



FIG. 10: The extra one loop term for the power spectrum with a non-Gaussian field. 
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FIG. 11: The two loop terms for the power spectrum which involve a 4-point function of the fields. 
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FIG. 12: Tree level terms of the 4-point function. 



of which are non-zero for a Gaussian field. Diagrammatically the terms involving a 3-point or 4-point function of 
the fields are given in Fig. 111! which corresponds to the first six terms of 



2 loop 



1 



c 



(27r)6 



dVrf''?2 ( iA^AiVBCDT^''^^(qi,q2 -qi,ki -q2,k2) 



+ -NabNcdT 
+NabNcdeB 



ABCD 



(qi,ki - qi,q2,k2 - q2) 



BDE 



(|ki + qil, |q2 - qil, |k2 - qsQC"^^ (gi) + ^NabNcdeB''''^ (k, qi, |k2 + qiQC^^fe) 



+ ^iVAA^BCDBS^''^(qi, |qi + qs | , ^2)^^'' (fc) + ^NaNbcdeB'''''' {qi, |ki - qi | , fc)C^^ fe) 



jCDE(„ i„ I „ I „ \r^ABii„\ I AT AT T^ABC/„ n, „ 1 i„\r<DE 1 

ANANBCDEFC^''{k)C''^{q,)C''^{q2) + ] Nabc NdefC^"" {k)C''^ iqi)C^'' iq2) 



+ lNABCDNEFC^''iqi)C^''{\\^i - q2|)C^^(g2) 

+ lNABcNDEFC^''{qi)C''^{\ci2 ^ qi|)C^^(|q2 - kil) 
6 



(35) 



The final four terms were already shown diagrammatically in Fig.lH they correspond to the last four terms of (j35p . 
However the mathematical expressions for these four diagrams given in (|16p are correct at leading order in slow 
roll, while the expression above is valid to all orders in slow roll. Note that 1 of the 1 loop diagrams and 4 of the 
2 loop diagrams have dressed vertices, these can be removed by renormalising the vertices in a similar way to the 
case with Gaussian fields, see App. lA 21 



B. Trispectrum 



The connected part of the primordial trispectrum is defined by 

(Cki Ck, Ck3 Ck4>c = rc(ki, k2, k3, k4)(27r)3<53(ki + k2 + ka + k4) 



(36) 



Diagrammatically all of the tree level terms for the 4-point function are given by Fig. 1121 Writing down these terms 
in order we get 



rf'=<=(ki,k2,k3,k4 



NaNbNcNdT^''^'' (ku-k^, k3, k4) 



+NabNcNdNe (C^^(fci)S^^^(lki -Hk2|,fc3,A:4) + (H perms)) 
+NabNcdNeNf (C"4^(|ki +k3|)C7^^(A:3)C^^(fc4) + (11 perms)) 
+NabcNdNeNf (C^^(fc2)C^^(A:3)C'^^(fc4) + (3 perms)) . 



(37) 



This result was first derived in [ij, |4ij without using diagrams (see also [231). Note that all of the numerical 
coefficients are one, because all of the terms are tree level, and that only the final two terms are non-zero if the field 
fluctuations are Gaussian. The real space diagrams for these two terms were given in Fig. [T] 
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VI. CONCLUSION 

We have presented a diagrammatic method for calculating any n-point function of the primordial curvature 
perturbation, at tree-level or any required loop level. Rules are given for first drawing diagrams at the required 
order, and then for writing down the corresponding terms in the rt-point function, in either real or Fourier space. 
For example, we have drawn diagrams corresponding to all the 1- and 2-loop terms to the power spectrum and have 
given the corresponding corrections to the power spectrum. 

Our method is based on using the (^A'^-formalism which identifies the local curvature perturbation with the 
perturbation of the integrated expansion from an initial hypersurface. The expansion on large scales (typically 
larger than the Hubble scale) can be determined as a (non-linear) function of the initial local values of scalar fields 
during inflation using the homogeneous (FRW) equations of motion in the separate universe approach. In practice 
we are working with fields that are smoothed on scales k > fcmax = (ai?)*. In particular we neglect any ultraviolet 
divergence of the fields on small scales. 

We show that it is possible to renormalise vertices by absorbing an infinite sum of terms, corresponding to loop 
corrections to vertices. These are automatically accounted for if we work in terms of the average local expansion 
{N{(f>{x))) and its derivatives, defined in Eq. for Gaussian fields, instead of the background expansion, iV((/)o). 
This removes terms such as {tp'^) which is divergent for scale- invariant spectra, leaving all n-point functions finite 
in real space so long as (iV(0(x))) and its derivatives are finite. However divergent terms remain in Fourier space 
due to integrals over loop momenta between vertices which may diverge in the infrared. A discussion of the effect 
of the large scale cut off on observables is given in [27l.[4^.^ 

Results for the primordial n-point function are given in terms of derivatives of N with respect to the fields at 
Hubble exit and the n-point functions of the fields evaluated at Hubble exit. 

In single- field models of inflation the curvature perturbation is constant on super Hubble scales, so it is possible to 
calculate the derivatives of N explicitly at Hubble exit during inflation [IJ] without needing to know the details of 
the end of inflation or reheating, etc. In this case both the bispectrum and trispectrum have been shown to be small, 
regardless of whether the fields are assumed to be Gaussian or not, but assuming the field has a standard kinetic 
term in the action. During inflation the derivatives of A'^ can be calculated explicitly even in multiple field models 
for a separable potential, this been done to second [3, IZsl . l46ll and third order [4l[ . Several specific multiple field 
models have been studied in greater detail, e.g. A/'-flation [47l l48j and a model of multiple field inflation in which 
A^ can be calculated exactly, [4^. For all examples considered so far in the slow roll regime the non-Gaussianity 
has proved to be small (see also [sol [sTI . [53|). however there are many models which explicitly generate a large 
non-Gaussianity after inflation. Examples include the curvaton scenario ■()] and modulated reheating (53| . also 
preheating, for example with hybrid inflation [53 |. In general a numerical calculation of the derivatives of N is 
required, the advantage of the 6N formalism is that we only need to solve the homogenous background equations. 

If one assumes the fields are Gaussian then only the 2-point function of the fields is required, but in general one 
requires the (r — l)-point function of the fields when working to r'th order for any n-point function of ^. So far 
the 3- and 4-point functions have been calculated, both at leading order in slow roll and ignoring an y q uantum 
mechanical loop corrections. It is possible but not expected that these corrections may be important [25l.l55j. In the 
most commonly considered case of slow-roll inflation, the field is extremely close to Gaussian at Hubble exit so it 
is a good approximation to assume that the field is Gaussian and that any observable non-Gaussianity is generated 
after Hubble exit. 

Methods to generate a large non-Gaussianity during inflation include having a break in the potential, [i^, or to 
have a non-standard kinetic term so that the non-Gaussianity of the fields is large at Hubble exit [5f|, for example 
DEI infiation inspired from brane world models, 
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Two very recent papers which discuss this in depth are j37l | by Lyth who estimates bounds on the size of the loop corrections to the 
2-, 3- and 4-point functions and 26, 43] by Seery who more explicitly calculates the size of the one-loop contribution to the power 
spectrum. He finds that any large loop correction will come from the loop correction to the curvature perturbation on super Hubble 
scales rather then the quantum mechanical loop corrections on sub Hubble scales. 
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APPENDIX A: RENORMALISING DERIVATIVES OF TV 

Here we provide a proof that renormalising vertices removes the diagrams with dressed vertices. In real space this 
is equivalent to removing terms with {</?^). We first provide a proof for the power spectrum where the notation is 
simpler. We then provide a proof for the general 71-point function, initially assuming that the fields are Gaussian at 
Hubble exit, Sec. lA 11 and finally extend this proof to allow for arbitrary field correlations at Hubble exit. Sec. lA 21 

The gauge invariant curvature perturbation is defined by ([2]) as 

where the subscript p means that there are p field indices above the brace. From the definition of the 2-point 
function is given by 

00 ^ 

p. (7—1 ^ ^ ' ^ ^ ^ 

P 1 

Because the fields are Gaussian every n-point function of the fields can be split into products of 2-point functions. 
We can choose the pairs to be connected pairs, i.e. each Lp is evaluated at a different point or disconnected pairs 
which are (1^9^}. We consider every possible way of splitting the term {(fi'^f^^ ■ ■ ■ fx'Py ' ' ' Vy) i^ito pairs: 

- Only disconnected pairs: this term vanishes from (jA2j because we are working with ^ which satisfies {(,) — 0. 

- 1 connected pair: there are pq different ways to select this term. This gives 

E ^,^AB ■■■C^DE--- fS'"'G{\x -y\){^^--- (Na) {N^)Gi\x - y\) , (A3) 

p q 

OC - 

where (iV^) = ^ ^__^iV^^ . . . ^(^^ . . . ^c:^ . (A4) 
^ ^ p-i 

- 2 connected pairs: there are p{p — l)q{q — l)/2! ways to choose these two connected pairs. The factor of 2! 
comes because the order in which we pick the pairs does not matter. This gives 



h ^ ip ^ 2)!(g - 2y. ''^B^N^E^s^''s''''G'i\^ - y\){- ■ ■ ■ ■ O 

p-2 q-2 

\NAB){N'''')GH\x~y\). (A5) 



2! 

Continuing this for every possible number of connected pairs, we find that (|A2[1 is equivalent to 



OC 

iCCy) = Y.-(^AB ■ ■ ■ c)iN^^-'')G''i\^ - y\) ■ (A6) 



r— 1 



1. General n-point function 

From the definition of C the n-point function is given by 

iUv-Q - E „,J...^, Nab---cNde---f---Ngh.--i 



nlol ■ ■ ■ r\ 
p,q,...,rt"^ 



X ((^^^f ■ • ■ - (^^^f • ■ • ^?)) • • ■ (^?^f ■■■vi- (^?^f ■ • • vi))) ■ (A7) 

Similar to the calculation of the 2-point function, we have to split the term {'P'^ip^ ■ ■ ■ ■ ■ ■ f'^^^ ■ ■ ■ ^i) into 
pairs in every possible way. In the case that some of the x,y, ■ ■ ■ z do not appear in a connected pair there is no 
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contribution to the n-point function, for example if ip^ terms do not appear in any connected terms then the last 
bracket of is ■ ■ ■ ^i) - {^^Vz ' ' ' ^i)) = 0- We are left with 

(CCy-C.) - E ^J^r—r,NAB---CNDE---F:--NGH---l 




x(^^^f •■•^^■••^^(pf ••.^^),,, (A8) 

where the subscript cc to the angular brackets means that they must be decomposed into pairs such that every 
x,y, . . . , z appears in at least one connected pair. We introduce the notation that Ixy > is the number of connected 
pairs of x and y and similarly for every other possible connected pair. Furthermore is the number of (^^j that 
form part of a connected pair, so = Ixy + ■ ■ ■ + Ixz- Note that we require 1 < Ux < p. 

We consider a term in (jASp from an arbitrary choice of l^y, ■ ■ ■ ,lxz- The number of ways to choose the Ixy 
connected pairs of x, y correlators is given by p{p — ^) ■ ■ ■ {p — Ixy + l)? ■ ■ ■ {q — Ixy + ^)/^xy^-, then the number of 
ways to choose the Ixz connected correlators linking x and z is given by (p — lxy){p — Ixy — 1) ■ • • (p — Ixy — Ixz + 
1) • • • r{r — 1) ■ ■ ■ {r — Ixz + ^)/lxz- Continuing this counting for all the pairs, we find the number of ways to split 
into the required number of connected pairs is 



p{p - I) ■ ■ ■ {p ~ Tlx + 1) ■ ■ ■ r{r - 1) ■ ■ ■ {r ~ Uz + I) 



1 1/ I . . . / I 

^xy-'^xz- ^yz- 



(A9) 



This choice of /^i/, • • ■ , lyz therefore gives the following term to (|A7I 



y. 7 ^^ 1 r,NAB-cNDE-F ■ ■ ■ Ngh-iS^'' 



xG'-(|x - 2/1) • • ■ G'- (|y - • • ■ (^) (AlO) 

{Nab . . . g) • • • i^GH ■ ■ ■ 1)5^'' 5^"" ■ ■ ■ G'- (|x - y\) ■ ■ ■ G'- (|y - z\) , (All) 



''xy- ''yz- 



where there are Ixy contraction of field indices of the first two factors of {Nab---c) etc. Repeating this for all allowed 
choices of ^a;^, • • • , lyz we get 

(C.C.---C.)= E , / ; , (iV^j^...g)---(iVgg...^)^^^^^^...G'-'^(|:r-y|)...G'«-(|y-z|). (A12) 

^xy . ■ ■ ■ l^yz ■ ' 



Diagrammatically a term G'^"" — ?/|) corresponds to Ixy propagators which all have both ends attached to the same 
spatial points x and y. This explains the numerical factor of Ixy^-^ as given by rule 5 in real space (or equivalently 
rule 6 in Fourier space). 



2. Extension to non-Gaussian fields 



Here we extend the proof that renormalising vertices eliminates the diagrams with dressed vertices also holds for 
non-Gaussian fields. By doing this we also derive the numerical factors that are associated with every diagram. 
The proof is more complex than before because we must consider every possible correlation of an arbitrary number 
of fields, rather than reducing every n-point correlation to products of 2-point functions. We now have 

{NaB-c) = NaB-C + \naB-CDe{^°v'^) + J^NAB-CDEF{v^v'^ip^) 

+ J^NAB-CDEFG{'f°V^V^V^) + • • ■ ■ (A13) 

The correlators here includes the connected and disconnected parts. 

Again the n-point function of C is given by (|A7|) and we can reduce it to (|A8|1 in a similar way to the case 
with Gaussian fields. However we now have to decompose the term in angle brackets into connected correlators of 
arbitrary order in every possible way. We extend the notation that Ixy is the number of terms like {fxfy) more 
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general correlators like Ixxzz is the number of terms like {ip'^Lp^ip'^ip^)c. Every I must have at least 2 distinct Xi in 
the subscript. We define rix as the number of x's contained in the correlators and again we require that 1 < < p 
and similarly for n^, etc. We consider the term in (|A8p coming from a particular choice of l^y, ■ ■ ■ ,lx---x---z: which 
we denote by 



= Yji fx ■ ■ Jfx - ■ ■fzipipzYc'-""'"' (A14) 
p w 

and the product on the last line is over all chosen correlators of the fields and we have also dropped the field indices 
on this line. We can repeat the calculation for the number of ways to split the total correlator of (jASP into the 
appropriate connected pairs. The only new point is that whenever an Xi appears more then once in the same lx---x---z 
the order in which we pick those x's does not matter, so we have to divide by the factorial of the number of x's 
raised to the power lx---x---z- In total we get 

p{p - I) ■ ■ ■ {p - Tlx + I) ■ ■ ■r{r - I) ■ ■ ■ {r - + I) 

Y{lx...x-z\{u\---w\Y^-- ' ^ ' 

where the product in the denominator is over all Z's and the number of XiH repeated in each individual I. 
Therefore this choice of Vs gives the following contribution to (jA7p . 



rn \T~\ \\i 5^ 7 ^i 1 rrNAB - cNDE - F ■ ■ ■ Ngh - i 

{[lx...x-zKu^- ■ ■ ■ w\y---- ^ ^^^ (p - nx)\ • ■ • (r - n,y. 

xP^^-^"{---^^) ■■■{■■■ ^i) (A16) 



Repeating this for all allowed choices of the I's this gives 

^^^^^ • • • ^ ..5 ^ ^ Wx...x..J(J---n.iy--^ ^''^R^ • • • (^G^)^""" ■ 

Eq. (jAlSp explains the numerical factor of rule 7 in Sec. |Vl except for the case of dressed vertices. The numerical 
factor for dressed vertices requires a calculation of the number of ways to reduce the n-point function in (jA13p into 
connected pairs of the required order. If / correlators of order n each dress the same vertex then we need to know 
the number of ways to split the ^n-point function into I lots of connected n-point functions. The result is 

where the denominator of l\ arises because the order in which we pick the n-point functions does not matter. The 
factor of {ln) \ is cancelled by the numerical factor in (|A13|) . The generalisation to more complicated dressed vertices 
is straightforward to calculate, for example the number of ways to split an (ni^i -I- n2Z2)-point function into li lots 
of connected ni-point fimctions and I2 lots of connected n2-point fimctions is 

1 (/ini+/2n2)! ^^20) 



h\h\ (ni!yi(n2!)'^ ' 

In Fig. [T3] we give two examples of these rules, both of which are 4 loop terms of the dressed vertex of Na, the 
first diagram has / = 2 and n = 3 while the second diagram has li = 2, ni = 2, ^2 = 1 and n2 = 3. The associated 
mathematical expressions are 

^ Nabcdefg{^^^'^^''){v''^''v'^) . (A21) 



2!3!2 

^ -,NABCDEFGH{^''^^){ip''^''){^''^'^^") . (A22) 



2!2!23! 



18 



•A 



/I /I 
\ I ' \ \ I 



FIG. 13: Two terms of the renormalised vertex (Na)- 
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